Inverse Problems and Imaging 
Volume X, No. OX, 200X, X-XX 



Web SITE: http://www.aimSciences.org 



A DIRECT SAMPLING METHOD 
FOR INVERSE SCATTERING USING FAR-FIELD DATA 

JiNGZHI Ll 

Faculty of Science, South University of Science and Technology of China 
Shenzhen, 518055, P. R. China 

JUN ZOU 

Department of Mathematics, The Chinese University of Hong Kong 
Shatin, Hong Kong 

(Communicated by the associate editor name) 

Abstract. This work is concerned with a direct sampling method (DSM) 
for inverse acoustic scattering problems using far-field data. Using one or few 
incident waves, the DSM provides quite reasonable profiles of scatterers in time- 
harmonic inverse acoustic scattering without a priori knowledge of either the 
physical properties or the number of disconnected components of the scatterer. 
We shall first present a novel derivation of the DSM using far-field data, then 
carry out a systematic evaluation of the performances and distinctions of the 
DSM using both near-field and far-field data. A new interpretation from the 
physical perspective is provided based on some numerical observations. It is 
shown from a variety of numerical experiments that the method has several 
interesting and promising potentials: a) ability to identify not only medium 
scatterers, but also obstacles, and even cracks, using measurement data from 
one or few incident directions, b) robustness with respect to large noise, and 
c) computational efficiency with only inner products involved. 

msg (ams 2010): 35r30, 41a27, 78a46 

1. Introduction 

The capability of effectively retrieving the location and/or geometrical features 
of unknown scatterers from the knowledge of the scattered wave field (near-field or 
far-field data) is of paramount importance in many practical applications such as 
underground mine detection, geophysical exploration in oil industry, detection of 
defects or cracks in nondestructive testing, target detection using radar or sonar sys- 
tems and ultrasound imaging in biomedical equipments [6l [12] . Both locations and 
geometrical shapes of scatterers are often the final target in the applications. Many 
quantitative methods have been developed for such purposes [3l [4[ [121 ITqI [22| [35 | [36] . 
But there is a crucial step for these imaging processes, to find some reasonable 
estimated locations and shapes as their initial computational sampling domains, 
otherwise the imaging processes may not work or work with huge computational 
efforts. 
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A promising direct sampling method (DSM) was proposed recently in [21 to 
recover the inhomogeneous medium scatterers from near-field measurement data 
for inverse acoustic medium scattering problems. Compared with other sampling- 
type schemes (see, e.g., [H [24l [32] for the detailed surveys) such as multiple signal 
classification (MUSIC) [i[Tl[l7l|20l|34] and linear sampling method (LSM) [aHOl 
[23] . a distinct feature of the DSM is its capability of depicting the profile of medium 
scatterers using one or very few incident wave fields, through the computationally 
very cheap inner products of the measured data and the fundamental solution in 
the homogeneous background medium. 

In this work, we will further explore in this direction by generalizing the DSM 
from the knowledge of near-field data to far-field data of the scattered waves for 
inverse medium and obstacle scattering problems. Extensive numerical experiments 
are done to investigate systematically the behavior of the DSM. An important new 
observation is that the DSM is able to recover not only medium scatterers, but also 
obstacles, and even cracks using one or two incident waves. 

The full inverse scattering problem (ISP) is to determine both the shapes and 
the physical profiles of the scatterers, hence may require the measured data at 
many frequencies, not the one from a single frequency as considered in this work. 
It is well known that the ISP is highly nonlinear, and the first key step for the 
numerical treatment of a nonlinear problem is to work out a good initial guess for 
its approximate solution. As DSM is computationally very cheap and works with the 
data from only one or few incident fields, it can naturally serve as a fast, simple and 
effective alternative to existing numerical tools for locating a reliable approximate 
position of the unknown scatterers, which can then serve as a good initial guess 
in any existing method (see, e.g., [H [12] [T9[ |22]) for achieving a more accurate 
estimate of the scatterer support and the inhomogeneity distribution. As we shall 
see numerically, the DSM can provide a very reliable location of each individual 
scatterer component so that one may start with a much smaller sampling region in 
a more accurate but computationally more demanding method. The reduction of the 
sizes of the initial sampling regions for unknown scatterers may save us an essential 
fraction of the computational efforts in the entire reconstruction process by most 
existing methods. Other approaches are available in the literature, which do not 
require an a priori good initial guess of the scatterers, e.g., the globally convergent 
numerical method. We refer to the monograph [5 and many references therein for 
the detailed description of the method and its theory for ISPs corresponding to 
some hyperbolic equations. Because of the time dependence, the ISPs of [5 can be 
considered as the ones with the data given at many frequencies. 

Inverse scattering reconstruction methods using one single incident wave, or sim- 
ply one-shot methods, date back to a long-standing open problem in the inverse 
scattering community, namely whether one can determine the unknown scatterer 
by using the near-field/far-field data from only one or several incident waves. The- 
oretically, the results on unique identifiability have been understood only partially, 
e.g., for some special classes of scatterers [12J[T3], or for the class of polygonal or 
polyhedral scatterers [2l [151 EH [SSl ISO] . Armed with the aforementioned unique- 
ness results, one-shot methods have been studied widely in the past few years, e.g., 
in [33l [Ml [211 HH] • This work extends the one in ^2]^ to make a more systemati- 
cal investigation of such a promising one-shot method applied to inverse scattering 
problems using near-field or far-field data for many possible cases of scatterers: ob- 
stacles, inhomogeneous media and cracks or their combinations. This study reveals 
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more potentials of the one-shot method and provides some physical hints to answer 
the pending open problem from the numerical perspective. We emphasize that the 
direct sampling method in the current work is based on the same indicator function 
as that in the orthogonality sampling method developed and studied in [33l [16] , 
but we shall present a different approach and motivation to derive this method and 
justify its effectiveness, in particular removing the "smallness" assumption that is 
crucial to the derivation of the method in [33^, . 

The paper is organized as follows. In Section [2] a brief review of inverse scattering 
problems is presented, along with some useful notations and identities. Section [3] 
describes the mathematical motivation of the DSM using far-field data and pro- 
poses a new indicator function. Section [4] provides extensive numerical simulations 
to evaluate the performance of the DSM using near-field or far-field data from ob- 
stacles, media and cracks. In addition, we shall provide a physical interpretation of 
the DSM, and compare the major features of the DSMs using near-field and far-field 
data. Some concluding remarks will be given in Section [5] 

2. A BRIEF REVIEW OF INVERSE SCATTERING PROBLEMS 

In this section we shall briefly describe the time harmonic inverse obstacle and 
medium scattering problem using near-field or far- field measurements pJJ 112j . Con- 
sider a homogeneous background space {N = 2,3) that contains some scat- 
terer components such as obstacles or inhomogeneous media, or both, occupying a 
bounded domain D. Let u'^^^ = ex.p{ikx • d) be an incident plane wave, with the 
incident direction d G S^~^ and the wave number /c, and u = u^^^ + be the total 
field formed by the incident and scattered fields. Then the total field u induced by 
the obstacles satisfies the Helmholtz equation 

(1) Au^k^u = in R^\D, 

or induced by the inhomogeneous medium scatterers satisfies 

(2) Au^k'^n'^{x)u = in R^, 

where n{x) is the refractive index. To account for the absorbing medium, the 
refractive index can be modeled by the complex form 

(3) = ni{x) + m2{x). 

The models above describe not only time-harmonic acoustic wave propagation, but 
also electromagnetic wave propagation in either the transverse magnetic or trans- 
verse electric models [121 EI] • 

Remark 1. Note that the obstacle scatterer can be viewed as the limiting case of 
the medium scatterer with vanishing or singular material properties. For example, 
an acoustic sound-soft obstacle is a limiting; case of (H as n2 +oo; see, e.g., I25| 
Eq. (4.4)] or [26l Sect. 4]. One may also refer to [27 for the sound-hard limiting 
case for general dimensions. Therefore, the indicator functions derived in the sequel 
are applicable to both inverse obstacle and medium scattering problems, although 
our derivations are based only on the medium scattering case. 

In the rest of this section we introduce some basic notation and fundamental 
functions that will be needed in the subsequent discussions. First, we define a 
coefficient function, 77 = (n^ — l) /c^, which characterizes the inhomogeneity of the 
concerned media and is supported in the scatterer D cR^. Then we define function 
/ = r]u^ which is called the induced current by the inhomogeneous media. Let 
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G(x, y) be the fundamental solution to the Helmholtz equation in the homogeneous 
background, that can be represented by (cf. [11] [T2]) 

-H^,'\k\x-y\) for TV = 2; 
(4) G{x,y) = ^ I exp(iA:|x-^|) 



, 47r k\x — y\ 



for TV = 3 



where Hq'^ refers to the Hankel function of the first kind and zeroth-order. With 
the help of the asymptotic properties of the fundamental solution G [12 , Eqs. (2.14) 
and (3.63)], we have for = 2, 3 that 

(5) »^(^) = |^^(lr'-!)j2 {"°°(^) + as \X\^00, 

with X = x/\x\ G S^~^. The scattered near-field and the far-field u"^ in ^ has 
the following very convenient representations (cf. [12, Eqs., (8. 12) and (8.27)]): 



(6) u'{x) = / G{x,y)I{y)dy, 

Jd 

(7) u^{x) = [ G^{x,y)I{y)dy, 



Id 

where the far-field pattern associated with the fundamental solution G is given by 

' exp(i7r/4) 



(8) G^{x,y) = { 



/Skn 



exp(— i/cf • y), N = 2; 



-^exp{-ikx-y), TV 



3. Direct sampling method using far-field data 

In this section we shall derive the DSM using far-field data. Like many other 
sampling- type methods such as LSM, MUSIC and factorization methods [T0^,^,'34], 
the essence of the DSM is to construct an indicator function which has significantly 
different behaviors inside and outside the scatterers, e.g., the indicator function of 
the LSM blows up outside the scatterer but remains finite within the scatterer. Let 
F be the surface where the near-field data is measured, 1] be a domain contained 
inside F such that the scatterer D lies in Q. For any sampling point Xp ^ the 
indicator function of the DSM using near-field data is given by (cf. [21 ) 

{u^ ^ G(-, Xp)) j^2(r) 

For the motivation and derivations of ([9| using near-field data we refer the reader 
to [21]. 

We are now going to derive an indicator function for the far-field data, a coun- 
terpart of ([9| for the near-field data. To do so, we first derive the following key 
lemma about the L^-correlation measure of the far-field data of two monopoles on 
the unit sphere. 

Lemma 3.1. For the far-field patterns G^ix^y) associated with the fundamental 
solutions G(x, y) (see (|8| and Q ), the following correlation holds for any two points 
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and the sampling domain Q: 

(10) / G^{x,Xj)G^{x,Xp)ds{x) = C^{G{xp,Xj)) 

where C is a constant depending only on the wave number k and the dimension N. 
Proof. By the definition of G^, it is easy to verify that 

(11) / G'^{x,Xj)G^{x,Xp)ds{x) = Ci / exp{ikx ■ {xp- Xj))ds{x), 

where Ci is a constant depending only on the wave number k and the dimension 
N. 

We carry out the proof in two steps, one for A/" = 2, and the other for N = 3. 

Step 1. For N = 2, (10) is a special case of Graf's addition theorem (TJ 
Eq. (9.1.79)] by showing an integral representation of Jn{kr) ex.-p{m6) for an in- 
teger n. To see this, we can write the right-hand side of (11) in polar coordinates 



(12) / exp(iA:x • {xp — xj)) ds{x) = / exp I ik{xp 

i§l J-TT \ 



Xj) 



cosO 
sin^ 



dO. 



It suffices to show that 
(13) Jn{kr) ex.-p{m6) 



/ exp ( ikx 



27r 



cos ( 
sine 



exp(in^ 



with r = and x = (rcos^, rsin^). Therefore we can arrive at (10) in two 



dimensions, by setting n = 0, replacing xhy Xp — xj in ( |13[ ) and using the fact that 

C^{G{xp,Xj)) = ^ (^iH^^\k\xp - Xjl)^ = Jo{k\xp - Xj\) . 
The formula (p^Sl) expresses regular cylindrical wave functions as a superposition of 



plane waves. To see (13), we first recall the Jacobi expansion [37^ Eq. (22)] 

oo 

(14) exp(i/crcos?/^) = ^ i^Jm{kr)ex.p{irmjj). 



m= — oo 



By replacing by 6> — in (14), then multiplying it by exp {—m{6 — (j))) and inte- 
grating over (— TT, tt) with respect to (j) we obtain 

(15) i^Jn(A:r) = ^ / exp (iA:r cos(l9 - ^)) exp (-in(l9 - 0)) d^, 

where we have employed the simple fact that exp(in?/^) dil) = 27r for n = and 
forn ^ 0. 

Now we can rewrite (p!5|) as 

(16) i^Jn{kr) 



exp {—'mO) 

2lT 



/' 

J — TT 



exp {ikr cos (j) cos 6> + r sin sin 0)) exp(in^) d^. 



which reproduces ( |13[ ) by a simple rearrangement. 

Step 2. For = 3, we use the following special variant of the Funk-Hecke 
formula (cf. [12, Eq. 2.44] and [31, p. 29]) 

(17) / exp(-i/cx • z)Yn{z) ds{z) = ^jn{k\x\)Yn{x), xeR\r>0 

^S2 1 

for spherical harmonics Yn{-) of order n and spherical Bessel functions jn{') of order 
n. Now we set n = and replace x hy Xp — Xj in (17). Being aware of the lowest 
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order spherical harmonics Yq{x) = l/v47r, which is a constant and thus can be 
dropped from both sides of the equation. Since jo{t) = sin(t)/t, we obtain that 

\ K\Xp Xj\ J 

in three dimensions. This completes the proof. □ 

Next, we divide the sampling domain Vt enclosing the scatterer D into a set of 
small elements {tj}, i.e., squares in 2D or cubes in 3D. Then we approximate the 
integral relation ^\ by the rectangular quadrature rule on each element: 

(18) = / G°°(i, y)I{y) dy « V u;,G°°(i, y,), 

Jo 

where the summation is over all elements Tj which intersect the sampling domain Vt^ 
the weight Wj is given by |Tj|/j, with |rj| being the area/ volume of the j-th element 
Tj [N — 2, 3) and /j the evaluation of I[x) at the center of Ty 



Now multiplying (18) by G^(x,Xp) for any sampling point Xp G Vt^ then inte- 
grating over the unit sphere S^~^ and using ([To|, we obtain 

(19) 



/ u'^{x)G'^{x,Xp)As{x) ^ Cy"wj^{G{xp,Xj)). 



As we illustrate later (see Fig. [6|, the right-hand side term ^(G(xp,Xj)) above 
approaches a constant when Xp tends close to some point scatterer xj and decays 
quickly as Xp moves away from xj. This behavior motivates us with the following 
important indicator function for any sampling point Xp G Q: 



(20) ^^(%) = n ooli 

where the L^-inner product occurring in the numerator is the integral given by the 



left-hand side of (19). In view of the singular behaviors of the fundamental solution 
G{x,y) in (19), the indicator ^"^{xp) approaches the unity when Xp tends close to 
or lies within the medium scatterer and decays quickly as Xp moves away from the 
scatterer; see Fig. |6] This indicating behavior will serve as the key ingredient in 
our new DSM using far-field data. 



We emphasize that the indicator function <l>^ in (20) was studied earlier in 
the orthogonality sampling method [33l [16] , but derived from a different approach 
and motivation, and under a "smallness" assumption. As it has been seen, our 
theoretical foundation does not require such smallness restriction, thus makes it 
possible to explain partially why the DSM works also well for thin or long scatterers 
like rings and cracks. This is confirmed by numerical tests (see Examples 6 and 7 
in Section [4|. 

Finally we make a short remark about computing the indicator function ^^{xp) 



in (20). Clearly the computing is purely explicit and direct, involving only a cheap 
scalar product and a normalization operation, and it does not involve any matrix 
inversion and any solution process. Moreover, normalizations are quite cheap too. 



In fact we can easily see that the denominator in ( 20 ) is actually a constant and does 
not depend on the sampling point Xp, so it can be computed once for all sampling 
points. This is quite different from the indicator function ([9| using near-field data. 
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Therefore, the DSM using far-field data is computationahy much cheaper than its 
counterpart using near-field data. 

4. Numerical examples 

In this section, we will carry out a systematic evaluation of the performance of 
the DSM using the far-field data for some inverse scattering benchmark problems, 
including scatterers of obstacle, medium and crack types. In particular, we shall 
also present the numerical reconstructions by the DSM using the near-field data for 
most examples so that we may clearly see the comparisons and distinctions between 
these two different types of measurement data. 

Let us first introduce our experiment settings. We shall present several examples 
to illustrate the applicability and robustness of the proposed method for determining 
the scatterers from both exact and noisy data. We shall take the unitary wavelength 
A = 1, and the wave number k = 27r. When it is not mentioned, only one incident 
direction d = (1, l)^/\/2 is employed in the examples; otherwise two incidents, 
di = (1, 1)^/^/2 and d2 = (1, — 1)-^/a/2, are used and mentioned explicitly. In all 
the examples, the scattered near-field data is measured at 50 points uniformly 
distributed on a circle of radius 4A centered at the origin, and the far-field pattern 
u'^ is observed from 50 uniform distributed angles on the unit circle S^. The noisy 
data and uf" are generated point-wise by the formulae 

(21) = u%x) + eCmax|^i^(x)| and uf = u'^ix) + e Cmax|^x^(x) | 

X X 

for both near-field and far-field data respectively, where e refers to the relative 
noise level, and both real and imaginary parts of the noise follow the standard 
normal distribution. Our near-field data are synthesized using the quadratic finite 
element discretization in the domain (—6, 6)^ enclosed by a PML layer of width 1 to 
damp the refiection. Local adaptive refinement scheme within the inhomogeneous 
scatterer is adopted to enhance the resolution of the scattered wave. The far-field 
data are generated approximately by the integral representation [12^, Eq. (3.64)] 
along the circle centered at the origin with radius 5 using the composite Simpson's 
rule: 

24 

5/(5, e)A9^^ Y.^hi + 4/2^+1 + /2i+2) 
15^ 

where /(r, 0) is the integrand of Eq. (3.64) in [12] in the polar form and fj is the value 
of / evaluated at the j-th quadrature node (r, 0) = (5, 2j7r/50), j = 0, 1, . . . , 50. To 
ensure the sufficient accuracy of the far-field data, we refine the mesh successively 
till the relative maximum error of the data from two successive meshes is below 
0.1%. The indicator function value is normalized between and 1, and the sampling 
domain is fixed to be 1] = [—2 A, 2A]^, which is divided into small squares of equal 
width = O.OIA. The contour plot of the indicator function value will be displayed 
as an estimate to the profiles of unknown scatterers. For multiple incident waves 
with directions (ii, (i2, • • • , we take the maximum of all indicator function values 
point- wise: 

(22) $^(xp) =max{^^(xp; di), ^^(x^; ^2), ••• ,^^{x^', d,)} , 

where ^^{xp] d^), / = 1, 2, . . . , i^, is the indicator function defined by ( [2Q| ) associated 
with the l-ih. incident direction di. We denote the DSM using near-field and far- 
field data by DSM(n) and DSM(f), respectively, for sake of comparisons of their 
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(d) (e) 

Figure 1. Example 1: (a) true scatterer; Reconstruction results 
using (b) exact near-field data, (c) noisy near-field data with e = 
20%, (d) exact far- field data, (e) noisy far-field data with e = 20%. 

performances. We shall carry out three groups of benchmark problems: medium 
scatterer, obstacle scatterer, and crack scatterer. 

DSM FOR MEDIUM SCATTERERS. We first test the four examples from [21 , but 
using both near-field and far-field data to compare the performances of two different 
types of data. 

Example 1 (A singular point medium). The example considers one square 
scatterer of width 0.02A located at the origin, with coefficient r] of the scatterer 
being 1. 

This point scatterer identification is a standard benchmark problem. The contour 
plots using the DSM(n) and DSM(f) are shown in Figure [l] With only one incident 
wave probing, the location of the point scatterer is accurately positioned even if the 
data is severely distorted. We observe similar behavior of both indicator functions 
using near-field and far-field data. Both DSM(n) and DSM(f) can locate the position 
of the "point" scatterer in an accurate and stable way. The presence of noise level 
up to 20% in both measurement data has little effects on the indicator functions. 

Example 2 (Well-separated medium scatterers) Two square scatterers of side 
length 0.3A, located at (— 0.8A, — 0.7A) and (0.3A, 0.8A) respectively, are considered, 
both with coefficient rj being 1. 

This example verifies that the DSM can capture multiple sources. The indicating 
contour plots of the DSM using near-field and far-field data are shown in Figure[2] 
Both DSM(n) and DSM(f) are able to identify the location of both scatterers even if 
the measured data is significantly perturbed. Through studying the contrast of the 
indicator contours of medium scatterers with respect to the surrounding background 
medium, one can choose the cut-off value ranging from 0.7 to 0.8 to truncate the 
profiles of the two scatterers. Choosing a large cut-off value helps reducing those 
misleading spurious medium scatterers, see yellow patches in Figure[2|b)-(e). 
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Figure 2. Example 2: (a) true scatterer; reconstruction results 
using (b) exact near-field data, (c) noisy near-field data with e = 
20%, (d) exact far- field data, (e) noisy far-field data with e = 20%. 



Example 3 (Two close medium scatterers). Here the same setting as Example 
2 is laid out, except that the two scatterers are now moved to (— 0.25A, 0) and 
(0.25A, 0), respectively. 

This example investigates the resolution limit when the DSM can separate mul- 
tiple sources. The results of near- field case and far-field case are shown in Figure[3] 
Both the DSM(n) and DSM(f) can separate these two close medium scatterers well. 
The location and size of the scatterers agrees well with the true ones even under 
large noise. When we further reduce the distance between those two scatterers (less 
than half a wave length), DSM cannot separate the two scatterer any longer, which 
is consistent with the Heisenberg uncertainty principle. 

Example 4 (A ring-shaped medium scatterer). We consider a ring-shaped 
square scatterer located at the origin, with the outer and inner side lengths being 
0.6A and 0.4A respectively and the coefficient 77 being 1. Two incident directions 
di = :^(1, 1)"^ and d2 = ~^)^ employed. The reconstruction results 

using near- field and far-field data are shown in Figures |4] and [5] respectively. 

This is a rather challenging task to estimate the profile of the scatterer. We see 
that the DSM(f) has similar performance to the DSM(n). One single incident wave 
is inadequate to revolve full features of the scatterer. Using two incident waves. The 
four corners of the square ring is correctly labeled by the peaks of the superimposed 
indicator function in Figure [4]^d), (g) and Figure|5|c), (f). 

From the above four examples, we see that the DSM(f) has similar performance 
as the DSM(n). Both methods work stably and can tolerate strong noise. These 
estimated locations and index values can serve as good initial guesses for further 
reconstructions of more accurate locations and inhomogeneity profile r]{x) of scat- 
terers through computationally much more demanding resolution processes such as 
nonlinear optimizations. 
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(d) (e) 



Figure 3. Example 3: (a) true scatterer; reconstruction results 
using (b) exact near-field data, (c) noisy near-field data with e = 
20%, (d) exact far- field data, (e) noisy far-field data with e = 20%. 



Physical interpretation. Before we study more challenging examples, we try 
to provide a physical justification of the DSM using both near-field and far- field 
data based on some numerical observations. 



Let us first consider the key factor ^ {G{xp^Xj)) in the correlation measure (19), 
which can be explicitly written (with an appropriate scaling factor) as: 



(23) C,v^(G(xp,x,)) 



{Mk\xj-Xp\), TV = 2; 
sin(/c|xj — Xp\) 



Jxi^X j 



A/- = 3, 



where Cat = 4 and 47r for = 2 and 3, respectively. Figure |6] shows that the leading 
term Cn ^ {G{xp^Xj)) of the correlation measure achieves the maximum when Xp 
approaches Xj, and decays quickly as the distance between them increases. 

Let us now first explain the DSM using far-field data. From Figure |6| it can 
be seen clearly that there is much stronger pattern correlation between u^{x) and 



G^(x,Xp) in (19) when r small, namly the sampling point Xp is 

sufficiently close to the center Xj of some inhomogeneous sampling element where 
Ij does not vanish due to the non- vanishing r]{xj). Let us take Example 1 to 
illustrate this key relation. We plot in Figure [7| the pointwise complex ratio of the 
far field data u^{x) to the far field pattern G^{x^ Xp) of the fundamental solutions 
located at the sampling point Xp being the origin. This complex ratio is computed 
with respect to 50 observation angles in [0, 27r] in the polar form with its polar radius 
and phase angle depicted in the left and right plots, respectively. It can be easily 
observed that both the polar radius and the phase angle in Figure [7| approaches 
some constants. In other words, u^{x) and G^{x^Xp) are almost linear dependent 
with each other. 
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(a) 






(e) 



(f) 



(g) 



Figure 4. Example 4: (a) true scatterer; reconstruction results 
using exact near-field data (b), (c), (d), and noisy data (e), (f), (g) 
with e = 20%; incident directions di, d2 and {di^d2) used respec- 
tively for (b) & (e), (c) & (f) and (d) & (g). 



For multi-component small-scaled medium scatterers with sufficiently large dis- 
tance between each other, multiple scattering is weak and thus ignored, so the exact 
scattered wave field can be approximated by its reduced counterpart, namely a 
superposition of small point sources induced by the incident plane wave illuminat- 
ing the scatterers. Thanks to the superposition principle of the wave phenomena, 
some scattered wave by an individual small medium scatterer shares common decay 
pattern with G{x^Xp) when Xp falls within the inhomogeneous medium, which leads 
to the parallelism of some component of u'^{x) and in the trace space 

on S^~^ as shown in Figure j?! Physically speaking, it amounts to exciting point 
sources when a plane wave impinges on some small-scaled media (size of media is less 
then a half wavelength). Then the far-field pattern G^{x^Xp) of the fundamental 
solution becomes a major component of the far-field u'^ when Xp is located within 
the scatterer. Therefore, if the indicator function approaches 1 for some 

sampling point Xp, it is plausible to claim that Xp is within the scatterer obstacle 
or medium. 

The reasoning for the DSM using near-field data can be viewed as a direct con- 
sequence of that of its far-field version. Again we use Example 1 for illustration 
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Figure 5. Example 4: reconstruction results using exact far- 
field data (a), (b), (c), and noisy data (d), (d), (f) with e = 20%; 
incident directions di, d2 and {di^d2) used respectively for (a) & 
(d), (b) & (e) and (c) & (f). 




(a) 



(b) 



Figure 6. Decay patterns of Cn ^ {G{xp^Xq)) as r = \xp — Xq\ 
increases in two dimensions (a) and three dimensions (b) when 
k = 27r. 



by taking the sampling point Xp to be the origin. Under the assumption of suffi- 
ciently large distance between the measurement surface F and the small scatterer 
inside (normally four or five times the wavelength), the traces of u^{x) and G(x, Xp) 
on F depend approximately linearly on their asymptotic amplitudes u'^{x) and 
G^(f,Xp), respectively, up to some complex scaling factors. This can be verified 
by checking the pointwise complex ratio of u^{x) to u'^{x) and that of G{x^Xp) 
to ^^{x^Xp), respectively, as showin in Figure [sj By recalling the nearly linear 
dependence of u'^{x) and G^(x, Xp) and the transitivity of the three linear approx- 
imations aforementioned, it leads naturally to the nearly linear dependence of the 
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Figure 7. Pointwise complex ratio ((a) polar radius; (b) phase 
angle) of u^{x) to G^{x^Xp) versus all the observation angle 0. 
Here Xp is chosen to be the origin, the center of the small inhomo- 
geneous scatterer in Example 1. 



traces of u^{x) and G{x^Xp) on F, which explain the effectiveness of the indicator 
function ([9| in [2T]. 

For multi-component scatterers, as in the far field case, the scattering field may 
well reduce to a superposition of small point sources induced by the incident plane 
wave illuminating the scatterers plus some high order terms due to the multiple 
scattering. When the sampling point Xp falls within a medium scatterer, one can 
obtain a nearly linear dependence of a significant component of the scattering field 

on the trace of the point source G{x^Xp) restricted on F. After normalization, 
such physical correlation of the near fields and G{x^Xp) leads to a fairly large 
value (nearly one) of the indicator function ^{xp) for the sampling point Xp. 



DSM FOR MIXED OBSTACLE AND MEDIUM SCATTERERS. The next example con- 
siders a mixed scatterer composed of an obstacle and an inhomogeneous medium 
with variable material coefficient. 

Example 5 (A mixed scatterer with obstacle and medium). The ex- 
ample considers two square scatterers of side length 0.3 A. The two scatterers are 
located at (— 0.8A, — 0.7A) (a sound-soft obstacle) and (0.3A, 0.8A) (an inhomoge- 
neous medium), respectively. The lower left one is an obstacle, which is excluded 
from the domain and denoted by a white square. We vary the coefficient r] in the 
upper right inhomogeneous square to study the wave interaction of obstacle and 
medium scatterers. 

This example verifies the capability of both types of DSMs to identify obstacles 
and also shows the condition under which the DSM can identify both obstacles and 
media. The results for = 1 and = 10 + lOi using the DSM(n) and DSM(f) are 
shown in Figures [9| and 10 respectively. For 77 = 1, it amounts to a refractive index 
n = \/l -\- l/(27r)^ for the right upper square medium scatter, which is very close 
to the background medium with unit refractive index. The scattering effect of the 
inhomogeneous medium scatterer is so small compared to the scattered wave due 
to the obstacle square. Therefore, one can only locate correctly the position of the 
obstacle and loses track of the information of the inhomogeneous medium scatterer. 
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2-K 



(a) 



(b) 



Complex ratio ((a) polar radius; (b) phase angle) of u^{x) to u^{x) versus the 

observation angle 0. 



(c) (d) 
Complex ratio ((c) polar radius; (d) phase angle) of G{x^Xp) to G^(x,Xp) versus 

the observation angle 6. 

Figure 8. Illustration of approximately linear dependence of the 
traces of u^{x) and G{x^Xp) on T on their asymptotic amplitudes 
u'^{x) and G^(f,Xp), respectively. Here Xp is chosen to be the 
origin, the center of the small inhomogeneous scatterer in Example 
1. 



However, when we increase the refractive index to be 10 + lOi, the multiple 
scattering interaction of the scattered wave field due to the obstacle and medium 
scatterers becomes more evident. Only in such circumstances the location of the 
medium scatterer emerges gradually as a red patch in the correct northeast location. 
In other words, the DSM performs well for mixed types of scatterers under the 
assumption that the interaction of scattered wave due to each individual scatterer 
(be it an obstacle or a medium scatterer) must be strong enough. Such good 
performance is quite robust with respect to large noise level. 

DSM FOR CRACKS. Now we are going to test two examples with cracks and demon- 
strate an interesting and promising application of the DSM to such an important 
scattering scenario. 
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(d) (e) 



Figure 9. Example 5 with 77 = 1: (a) true scatterer; reconstruc- 
tion results using (b) exact near-field data, (c) noisy near-field data 
with e = 20%, (d) exact far- field data, (e) noisy far-field data with 
e = 20%. 




(c) (d) 



Figure 10. Example 5 with = 10 + lOi, : reconstruction 
results using (a) exact near-field data, (b) noisy near-field data 
with e = 20%, (c) exact far- field data, (d) noisy far-field data with 
e = 20%. 
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(d) (e) 

Figure 11. Example 6: (a) true crack scatterer; reconstruction 
results using (b) exact near-field data, (c) noisy near-field data with 
e = 20%, (d) exact far-field data, (e) noisy far-field data with 
e = 20%. 



Example 6 (A horizontal crack). This example considers a thin crack of 
thickness 0.1 A, centered at the origin and parallel to the x-axis lying in the middle 
of the domain with length 1. The medium scatterer with coefficient 77 = 1 is illumi- 
nated by a horizontal incident plane wave, see Figures [iTJ a) for the configuration. 

Cracks are the most illusive scatterer type to be identified as it has a very small 
thickness, and it is non-trivial even with the data from multiple incident directions. 
The results using the near-field and far-field data from just one incident field are 
shown in Figure [TT] 

In the noise-free cases, both DSM(n) and DSM(f) can successfully determine the 
location and length for the crack, and the identified geometry is roughly a red long 
bar with correct length. It is important to observe that the DSM also works for 
this difficult example with one incident, even with large noise, say e = 20%; see 
FigurepT] (c), (e). But in this case, the DSM(f) performs better than the DSM(n), 
which seems to break the crack into two pieces. 

Example 7 (An L-shape crack). The example considers an L-shape crack 
with thickness 0.1 A and length 2A. The coefficient rj in the inhomogeneous region 
is 1. 

This L-shaped crack is more challenging than the crack in Example 6. We use 
the incident wave along direction d2 = (1, —1)^/V2 and the reconstruction is given 
in Figure 12 One can see that an L-shaped dark yellow bar is identified by choosing 
the cut-off value about 0.7 in the noise- free case for both DSM(n) and DSM(f). But 
with 20% noise, the DSM(f) provides a much better profile than the DSM(n) in 
terms of the length and connectedness of the crack. 

Finally we test with two incident wave directions di and and report the results 



in Figure 13 One can see that both DSM(n) and DSM(f) can yield significantly 



enhanced images than those obtained using a single wave in direction di or d2. 
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(d) (e) 



Figure 12. Example 7 with incident direction c^2 = (l, -1)^/V^ 
: (a) true L-shape crack scatterer; reconstruction results using (b) 
exact near-field data, (c) noisy near- field data with e = 20%, (d) 
exact far- field data, (e) noisy far-field data with e = 20%. 




Moreover, the recovered L-shaped geometry degenerates rather slowly as we increase 
the noise level from 5% to 20%, see Figure 13 This example demonstrates the 



robustness and effectiveness of the DSM as a potential promising technique for 
crack detections. 
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5. Concluding remarks 

A systematic evaluation of the performance of the DSM using both far-field data 
and near-field data has been carried out for some inverse scattering benchmark 
problems, such as scatterers of obstacle, medium and crack types. The numerical 
simulations have demonstrated the robustness and effectiveness of the DSM to iden- 
tify the locations of obstacles, inhomogeneity media, as well as thin cracks, with just 
one or two incident directions. An important advantage of the DSM is its tolerance 
against large noise in the observation data. 

The DSM method involves only simple operations like inner products, without 
any inversion or solution process required, so it is computationally very cheap, hence 
not intended for an accurate reconstruction of unknown scatterers and their physi- 
cal inhomogeneity. As it is cheap and works with the data from even one incident 
field, it can serve as a fast, simple and effective numerical tool for locating reliable 
approximate positions of unknown scatterers. The reconstructed shapes of scatter- 
ers and the final numerical indicator function values by DSM may then serve as a 
good initial guess for any other more advanced method to achieve a more accurate 
estimate of the scatterer shapes and the inhomogeneity distribution. In addition, 
with the reliable location of each individual scatterer component extracted by the 
DSM, one may start with a much smaller sampling region in a more accurate but 
computationally much more demanding method (see, e.g., [H |12[ [T9[ [22]). Con- 
sidering the severe ill-posedness of the inverse scattering problems, the reduction 
of the sizes of the initial sampling regions for unknown scatterers may save us an 
essential fraction of the computational efforts in the entire reconstruction process. 
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